ON THE EXTREMAL RAYS OF THE CONE OF POSITIVE, POSITIVE 

DEFINITE FUNCTIONS 



PHILIPPE JAMING, MATE MATOLCSI & SZILARD GY. REVESZ 

Abstract. The aim of this paper is to investigate the cone of non-negative, radial, positive- 
definite functions in the set of continuous functions on R d . Elements of this cone admit a 
Choquet integral representation in terms of the extremals. The main feature of this article 
is to characterize some large classes of such extremals. In particular, we show that there 
many other extremals than the gaussians, thus disproving a conjecture of G. Choquet and 
that no reasonable conjecture can be made on the full set of extremals. 

The last feature of this article is to show that many characterizations of positive definite 
functions available in the literature are actually particular cases of the Choquet integral 
representations we obtain. 
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1. Introduction 

Positive definite functions appear in many areas of mathematics, ranging from number 
theory to statistical applications. Since Bochner's work, these functions are known to be 
characterized as having a non-negative Fourier transform. 

Before going on, let us first fix some notations. We will define the Fourier transform of a 
function / G L 1 ^) by 

= /(£)= / f(x)e 2l ^dx 

and extend this definition both to bounded measures on R d and to L 2 (R d ) in the usual way. 
Here (•, •) is the scalar product on R d and | • | the Euclidean norm. 

A continuous function / is said to be positive definite if, for every integer n, for all 
x\, . . . ,x n G M. d , the n x n matrix [f(xj — Xk)]i<j,k<n is positive definite, that is, if 

n 

(1) y~] CjCj^f(xj -x k )>0 for all ci,...,c n G C. 

j,k=i 

Then Bochner's Theorem [Boc| shows that / is positive definite if and only if / = fl for 
some positive bounded Radon measure on W d (a probability measure if we further impose 
/(0) = 1). There are many proofs of Bochner's theorem, the nearest to the subject of this 
paper being based on the Choquet Representation Theorem, due to Bucy and Maltese [BMJ , 
see also \Be\ IChll IPh] . Let us recall the main features of these proofs, and thereby also some 
definitions. An element / of a cone Q B f is an extremal ray generator of O (or simply an 
extremal) it f = fi + fa with fa, fa E fl implies fa = A/, fa = fj,f, A,/i > 0. The first step 
in the proof of Bochner's theorem is then to show that the characters e 2i ^ x '^' are the only 
extremal rays of the cone of positive definite functions. The second step is to show that the 
cone of positive definite functions on R d is well capped, i.e. is the union of caps (compact, 
convex subsets C of f2 such that 0, \ C is still convex). It then follows from the work of 
Choquet that every element of such a cone is an integral over extreme points with respect to 
a conical measure. For more details, we refer to the references given previously. 

Let us now note that Bochner's Theorem, though being powerful when one wants to 
construct positive definite functions, may be difficult to use in practice. This mainly comes 
from the fact that explicit computations of Fourier transforms are generally impossible. For 
instance, it is not known precisely for which values of A and k the function (1 — |x|^)i is 
positive definite on M. d . This problem is known as the Kuttner-Golubov problem and we refer 
to |Gn2] for more details and the best known results to date. To overcome this difficulty, one 
seeks concrete and easily checkable criteria that guarantee that a function is positive definite. 
The most famous such criterion is due to Polya which shows that a bounded continuous 
function on R which is convex on [0, +00), is positive definite. More evolved criteria may be 
found in the literature (see Section [5] for more details). 

As it turns out, the functions so characterized are not only positive definite but also non- 
negative. We will call such functions positive positive definite. Such functions appear in 
many contexts. To give a few examples where the reader may find further references, let us 
mention various fields such as approximation theory [BuJ, spatial statistics [Bon], geometry 
of Banach spaces [Ko] , and physics |GS| . Despite a call to study such functions by P. Levy in 
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|Lej they seem not to have attracted much attention so far. To our knowldge, there are only 
two papers specifically devoted to the subject in the litterature: [GP] which is motivated by 
applications in physics and the (unpublished) paper [Borj which is motivated by problems in 
number theory. 

Before going on with the description of the main features of this paper, we will need some 
further notations. For d > 1 we define C r (M. d ) to be the space of radial continuous functions, 
and we stress the fact that in the sequel we will only consider radial functions in dimensions 
higher than 1. Now, let 



topology). For / 6 we denote by 1(f) = STj~n(/ — STj~) the interval generated by /. Then 
/ is an extremal ray generator if /(/) = {A/ : < A < 1}. As the cone is well-capped, 
so is Qj. Therefore, Choquet Theory applies and every positive positive definite function 
admits an integral representation over extremals. It is therefore a natural task to determine 
the extremals of 

We are unfortunately unable to fulfill this task completely. One difficulty is that among 
the extremals of the cone of positive definite functions, only the trivial character 1 is still 
in the cone of positive positive definite functions and is of course an extremal of it. We are 
nevertheless able to describe large classes of extremals. One such class is included in the 
compactly supported extremals of Ojr. In this case, we show that if an element of with 
compact support has a Fourier transform whose holomorphic extension to C has only real 
zeroes, then this element is an extremal ray generator. One would be tempted to conjecture 
that this describes all compactly supported extremals, but we show that this is not the case. 
Nevertheless, this theorem allows to show that many examples of functions used in practice 
are extremals, as for instance the functions (1 — |x| 2 )+ * (1 — |x| 2 )" for suitable ct's. 

The next class of functions we investigate is that of Hermite functions, that is, functions of 
the form P(x)e~ Xx , P a polynomial. This is a natural class to investigate since the elements 
of the intervals they generate consist of functions of the same form. This will be shown as a 
simple consequence of Hardy's Uncertainty Principle. Further, a conjecture that oral tradition 
|Goj attributes to G. Choquet (although P. Levy may be another reasonable source of the 
conjecture) states that the only extremals are the gaussians. This conjecture is false, as our 
results on compactly supported extremals show. We will construct more counter-examples 
by describing precisely the positive positive definite functions of the form P(x)e~ nx where 
P is a polynomial of degree 4 and showing that this class contains extremal ray generators. 

2 

This further allows us to construct extremals of the form P(x)e~ 7TX with P polynomials of 
arbitrary high degree. 

Finally, we show that most (sufficient) characterizations of positive definite functions ac- 
tually characterize positive positive definite functions and are actually particular cases of 



n d = {/ G C r (R d ) : f is positive definite} 
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Choquet representations. More precisely, it is easy to see that if <p is an extremal ray gener- 
ator in then so is ft{x) = (p(tx). It follows that 



is a positive positive definite function (for suitable [i) that is obtained by a Choquet repre- 
sentation with a measure supported on the family of extremals {(ft}- For instance: 

- as a particular case of our theorem concerning compactly supported functions, we 
obtain that the function (p(x) = (1 — |x|)+ is extremal. Polya's criterium may be seen as a 
characterization of those functions which may be written in the form ([2]) . 

- a criterium for deciding which functions may be written in the form ([2]) with (p(x) = 
(1 — |x| 2 )-|- * (1 — |a;| 2 )+ has been obtained by Gneiting. 

The article is organized as follows. In the next section, we gather preliminaries on positive 
definite functions. We then turn to the case of compactly supported functions in Section O 
followed by the case of Hermite functions. We then devote Section [5] to link our results with 
various criteria available in the literature. We conclude the paper with some open questions. 



2.1. Basic facts. 

Fact 1 (Invariance by scaling). 

Let f be a continuous function and A > 0. We write f\(x) = f(Xx). Then f G ST^ if and 
only if f\ G fl^ ■ Moreover f is an extremal ray generator if and only if fx is extremal. 

Another consequence is the following lemma: 

Lemma 2.1. Let uj G and v be a •positive bounded measure on (0, +oo). Define 



Then F G Q^. Moreover, write uj{x) = ojq{\x\) and assume that ujq is non-increasing on 
(0, +oo) and (strictly) decreasing in a neighborhood ofO. Then F is an extremal ray generator 
if and only if uj is an extremal ray generator and v = 5 a is a Dirac mass for some a G (0, +oo) . 

Proof. Note that the fact that F is continuous follows from Lebesgue's Theorem since uj is 
bounded and continuous. It is then also obvious that F is positive and positive definite. 

If ui is an extremal ray generator and v = 5 a then F(x) = uj{x/a) is clearly an extremal 
ray generator. It is also immediate that if uj is not extremal then F is not extremal. 

Finally, assume that the support of v contains at least 2 points a < b = a + 3n and let 
ifj be a continuous non-increasing function such that < ip < 1 and ip{x) = 1 on (0, a + rj) 
while ^ = 0on (a + 2r], +oo). Define 



(2) 




2. Preliminaries 





The measures tp{t) &v{t) and (l - ip(t)) &v{t) are bounded, so that F\,F2 are continuous, Fi 
and i*2 are positive positive definite, F = F\ + F%. 
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Assume now towards a contradiction that F\ = XF with < A < 1. From the assump- 
tion on ip, we easily deduce that < A < 1. Let to be a solution of ifi(t) = X. From 

u(x/t)(ip(t) - A) du(t) = we get that 

rto r+oo 

(3) / Ld(x/t)(ip(t) - X) dv(t) = / w(x/t)(X-i/f(t))di/(t). 

JO Jt 

In particular, for x = 0, we obtain 

"to r+oo 

(m~ A) du(t) = / (X-m)dv(t). 

Jto 

Now, X(o,t ) yPit) ~ A) di/(t) and X(i ,+oo)(A — ^>(i)) di/(i) are positive non-zero measures. 
Then, for |xq| > small enough to have that u> is strictly decreasing on (0, |a?o |/*o) 5 we get 



° u(x /t)(^(t) - X)du{t) = f° u (\x \/t)(ijj(t) - X)du(t) 

Jo 

< Lo (\x \/t ) / °ty(i)-A)di/(t) 
Jo 

r+oo 

= M\xo\/t ) / (A - 4>{t)) di/(i) 

Jto 

r+oo 

< / wo(|o;o|/t)(A-V(t))di/(t) 

Jto 

-oo 

u{x /t){X-ip(t))du(t) 

to 

a contradiction, so that F\ ^ XF and F is not an extremal ray generator. □ 
Fact 2 (Invariance under products and convolution). 

If f>9 € th en fd ^ ■ Further, if f,g are also in L 2 (say) then f * g £ Slf. 

Further if either f or g (resp. f or g) does not vanish, then , for fg (resp. f * g) to be 
extremal, it is necessary that both f and g are extremal. 

Indeed, assume that g is not extremal and write g = gi + g2 with gi/g,g2/g not constant, 
then fg = fg\ + fg 2 . Now, if fgi = Xfg then gi = Xg on the support of /. 

The converse is unclear and probably false. A possible counter-example may be constructed 
as follows. Assume there is a compactly supported extremal / such that / 2 is also extremal. 
Without loss of generality, we may assume that / is supported in [—1,1]. Let g = (45o + 
5-4tt + 5-4-n-) * f then fg = 4/ 2 would be extremal but g is not. 



2.2. Bessel functions and Fourier transforms. 

Results in this section can be found in most books on Fourier analysis, for instance \Gr\ 
Appendix B]. 
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Let A be a real number with A > —1/2. We define the Bessel function J\ of order A on 
(0, +oo) by its Poisson representation formula 

Let us define J-i/2(x) = cosx and for A > —1/2, J\{x) := Jx ^ . Then J\ extends to an 
even entire function of order 1 and satisfies J\(x) real and J\(ix) > for all x £ M.. It is 
also known that J\ has only real simple zeroes. 

As is well-known, if / is a radial function given by f(x) = /o(|x|), then its Fourier transform 
is given by 

/(6 = ^_i(/o)(|£|) 

2 

where 

r+oo 

Mf )(t) = (2ir) x+1 / / (r) J A (2vrrt)r 2A+1 dr. 
J o 

Bochner's Theorem has been extended to radial continuous positive definite functions by 

Schoenberg [Scl page 816] (see also |SvP| ) : a function ip is radial positive definite and radial 

on W 1 , d > 2 if and only if there exists a positive bounded measure \x on (0, +oo) such that 

(p{x) = <po(\x\) with 

i f+°° d r+oo 

(po ( r )= r -2+ 1 Jd_ 1 (2vrrs)d/i(s) = (2tt)2- 1 / Ja ^irrs) 52 _1 dju(s). 

Jo 2 Jo 2 

For = 1 this coincides with Bochner's theorem. 

We will also use the following well-known result: denote = t or according to t > 
or not. Let m a (x) = (1 — |x| 2 )". Then 

^(^ r( ;: 1)j y' ) -^r (a+1 ) J{+ „ (2 ^, 

3. Compactly supported positive positive definite functions 

In this section, we consider compactly supported positive positive definite functions. It is 
natural to look for extremals inside this class of functions because of the following (trivial) 
lemma : 

Lemma 3.1. Let f be a continuous radial positive positive definite function with compact 
support. Then the interval 1(f) contains only positive positive definite functions with support 
included in supp/. 

Moreover, if we write B(0,a) for the smallest ball containing supp/ ; i.e. convsupp/ = 
B(0, a), and if f = g + h with g,h C 1(f), then at least one of conv suppg and conv supp h 
is B(0,a). 

Proof. First if g G /(/) then < g < f so that supp g C supp/. 

Next, as / is radial, there exists a such that convsupp/ = B(0,a). Assume now that 
/ = g + h with g,h £ /(/)• As g (resp. h) is radial, the convex hull of its support is a ball 
and we denote it by B(0,b) (resp. B(0,c)). As g and h are both non-negative, the convex 
hull of the support of g + h = f is then B(0, max(6, c)) = B(0, a) thus the claim. □ 
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We may now prove the following: 

Theorem 3.2. Let f G C(M. d ) be a compactly supported positive positive definite radial func- 
tion. Write f(x) = fo(\x\) where fo is a compactly supported function on M + so that J±_-> (fo) 

2 

extends analytically to C. 

Assume that J±_, (fo) has only real zeroes, then f is an extreme ray generator in the cone 

2 

of continuous positive positive definite radial functions. 

Moreover, assume that Jd_, (fo) has only a finite number N of non-real zeroes, and let 

2 1 

g G 1(f)- If we write g = go(\x\), then Jd_, (go) has at most N non-real zeroes. 

2 1 

Proof. Let us write / = g + h with g, h radial positive positive definite functions. Then g and 
h are also compactly supported. 

Write g(x) = go(\x\) and h(x) = ho(\x\). It follows that J&_* (fo), Jd, (go) and J±_, (ho) 

2 2 2 

all extend to entire functions of order 1. From Hadamard's factorization theorem, we may 
write Jd_A (fo) as 

2 

(4) Jd_ 1 (fo)(z) = z k e az+b [J i 1 -?)^? 

where 2(f) is the set of non-zero zeroes of l 7a_ 1 (fo)- Let us further note that 

2 

(i) Ji_ x (fo) is real if z G R, thus a, b G R and if C G Z(f) then £ E Z(f); 

(ii) J7di (fo) is non-negative if z G M, thus real zeroes are of even order; 

2 

(hi) J|_ x (/o)(0) =f(0) = J / / since / > thus k = 0; 

(iv) Jd , (fo) is even thus if C G ^(/) then — £ G and a = 0. 

2 1 

We may thus simplify © to J^C/bXz) = f(0)E f (z)P f (z) with 

w= n (i-S) 2 

where = Z(f) D (0, +oo) and 

w= n 

C€2q(/) 

where Zn(f) = {( £ ^(/) : ReC > &; ImC > 0}. Similar expressions hold for Jd_, (go) 

2 

and Jd_ 1 (ho). It should also be noticed that both Ef and Pf are non-negative on the real 
and the imaginary axes. 

Now, let us assume that i7d_ 1 (/o) has only finitely many non-real zeroes, so that Pf is a 

polynomial. In particular, there exists an integer N and a constant C such that < 
C(l + \z\) N . 

As Jd_ x (/ Q ) = Jd^So) + Jd_ x (h ) with Jd.^/io) > 0, we get < Jd_ 1 (g )(z) < 

2 2 2 2 2 

^J-ii/o)^) for ^ rea/. It follows that Z + (f) C Z + (g), with multiplicity. Thus, we may 

2 

partition the multiset Z + (g) = Z + (f) U Z'(g). 
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From Hadamard's factorization, it follows that 

Ji_ 1 (go)(z)=g(0)Ef(z) ]J (l-*Yp g ( z ). 

Ce2'{g) V * 7 

So, we have written Jd_ 1 {go){z) = G(z)Ef(z) where G is an entire function of order at most 
1. Note that < (g )(x) < Jd, (fo)(x), so that < \G(x)\ < C(l + 1x1)^ for x real. 

2 2 

Further, from the positivity of / and g, 

2-irtx a„ ^- I t(„\„2-Kt% 



(5) < J*_ x (g*)(it) = / g(x)e Mx dx < / f( X y™ dx = J^ifoW). 

2 JB(0,a) JB(0,a) 2 



It follows that |G(z)| is also bounded by C(l + \z\) N on the imaginary axis. By Phragmen- 
Lindelof's Principle G is bounded by 2 N C{1 + 1^1)^ over each of the four quadrants Q ei , e2 = 
{eiRez > &: £2^1 z > 0}, e\ = ±1,£2 = ±1- From Liouville's Theorem, we thus get that 
G is a polynomial of degree at most N. 

In particular, if N = 0, then G is a constant and Jd, (go) = AJd_, (/o) thus g = A/. It 

2 2 _ 

follows that h = (1 — A)/, and / is an extremal ray generator. □ 

Note that in the course of the proof we have shown that if the Fourier transform / of a 
compactly supported positive positive definite function / has a non real zero £, then £ is not 
purely imaginary (see (0)) and — (, £, — £ are also zeroes of /. 

It should be noted that checking whether a particular function is an extremal positive 
positive definite function may be difficult in practice. Neverteless, we will now give a few 
examples. 
Example. 

Let us consider the characteristic (indicator) function X[-i,l] °f the interval [—1,1] in R. 
Then p(x) := * X[— 1,1] i s an extremal ray generator since its Fourier transform is 

m = 

The celebrated positive definiteness criteria of Polya characterizes those functions that 
may be written in the form J p{x/r) dfi(r) with [i a positive measure (see Section [5.ip . 
Further examples are obtained by convolving p's: 

Xl-nAn/2] * X[-n/2,ri/2] * ■ ■ ■ * X[-r k /2,r k /2] * X[-r k /2,r k /2]- 

Another class of examples is given by the following: 
Corollary 3.3. For a > —1/2 define the function m a on M. d by 

((1 - \x\ 2 ) a if\x\ < 1 



2 



m a (x) = (1 — \x 



2\a 



otherwise 



Then, m a * m a is a continuous positive positive definite function that is an extremal ray 
generator. 

This covers the previous example and extends it to higher dimension since mo is the 
characteristic (indicator) function of the unit ball of M, d . 
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Remark. If we were considering positive definite tempered distibutions, then this result stays 
true for a > — 1. 

Proof. First of all, note that if a > 0, m a £ L°° and for — 1 < a < 0, m a £ L p for all 
a > —1/p. It follows that m a * m a is well defined, has compact support and is continuous 
for all a > -1/2. 

Further, as is well-known (see e.g. \Gi\ Appendix B5]): 

- (C\ - r(Q + l) Jd/2+a(27T|g|) 
m ^> ~ ^ a |£|d/2+ a ■ 

It follows that the Fourier transform of m a * m a is given by 

m _ r(a + i) 2 J d/2+a (Mi\? 

m a * m a {^) - ^ |^|d+ 2 Q • 

But, from a Theorem of Hurwitz [Hu] (see \Wa\ 15.27 page 483]), Jd/2+a nas only positive 
zeroes since d/2 + a > —1. The corollary thus follows from the previous theorem. □ 

Example. 

A particular case of the previous result is that the function w(x) := m\(x) * m\(x) = (1 — 
x 2 )+ * (1 — x 2 ) + on M is an extremal ray generator. The function w has been introduced in 
the study of positive definite functions by Wu [Wuj in the context of radial basis function 
interpolation. A simple criteria for writing a function in the form J w(x/r) d/i(r) with \jl a 
positive bounded measure has been recently obtained by Gneiting [Gnl] (see Section [5^2]) . 

Further examples are then obtained by taking scales w r of w and convolving several w^s 
together: 

(1 - x 2 /rl) + * (1 - x 2 /rl) + * ••• * (1 - x 2 /r 2 k ) + * (1 - x 2 /r 2 ) + 

Remark. 

When A — > oo we have m\(^/ : Kx/\f\) — ► e _7nr pointwise, in L 2 , and uniformly on compact 
sets. 

It would be tempting to conjecture that every compactly supported extreme ray generator 
is covered by Theorem 13.21 i.e. it is of the form that its Fourier transform has only real 
zeroes. Nevertheless, this is not the case: 

Proposition 3.4. There exists a continuous compactly supported extreme ray generator f of 
£lf such that f has non-real zeroes. 

Proof. Consider <p(x) = (1 — x 2 ) 2 ^ * (1 — x 2 ) 2 ^. A cumbersome computation (by computer) 
shows that 



ip(x) 



65o(2 - xf(x A + 10|x| 3 + 36x 2 + 40|x| + 16) if |x| < 2 
otherwise 



Let r > and < 9 < tt/2 and £ = ^-e ld , let us define f r e = f( by 
(6) = (l + ^^(2.vre) 2 + l(2 i vr0 4 )^). 
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Then is positive definite. 

Further, as p is smooth, and as d k p{£) = (2i7r£) fe </?(£) we get that 



, . , / 2cos20 o2 1 A , , 
f c (x)= ^l+_ 1 _^ + _ 9 4j ^ {x y 



The computations are easily justified by the fact that (p is of class C 4 . Note that, for |x| < 2, 



and 



(f/'(x) = — (3x 4 + 18|x| 3 + 30|x| 2 - 12|x| - 8)(2 - x) 3 
105 



ipW(x) = -(3x 4 + 6|x| 3 - 8x 2 - 16|x| + 8)(2 - x). 
5 



We will now take 9 = 7r/4. Actually, we believe that for every < 6 < ir/2 there is a 
unique r such that is extremal. 
We then have 

f ( (x)=p(x)+p^(x)/r 4 . 
It is not hard to see that <p is decreasing and positive on (0, 2) and that tpW is decreasing on 



[0,J2- 2/V3) and on (y 2 + 2/ \/3, 2) increasing on (y 2 — 2/ \/3, y 2 + 2/ \/3), positive on 
(0, xi) U (x2,2) and negative on (xi,X2) where x% = 0.441... and X2 = 1.462.... In particular, 
for r big enough (r = 4 will do), f r ^u is positive on (0, 2) and for r small enough (r = 3 will 
do), has two zeroes on (0,2). Therefore, there exists a unique r such that / rjVr /4 has 

exactly one double zero on (0,2). A numerical computation shows that r ~ 3.342775. Let us 
denote this zero by x^. A computer computation shows that x^ ~ 1.303. 

We will need a bit more information. Let £ = ^e*^. Assume that > and that there 
exists xg G (0,2) such that f^(x^) = 0. Then, as /c is a polynomial on (0,2), we also have 
f'Axg) = that is 



(7) 



f ^) + 2 ^^(x ? ) + ^^ 4 )(x c ) = 

\^(x 5 ) + 2^^ (3) (^) + > (5) (^) = 0- 

A computer plot shows that ip"(p( 5 ' — ip^'ip^ ^ on [0,2) so that, taking the appropriate 
linear combination of both equations, we get 

p CQS 2ip ^ ip'jx^ip^jx^ - p{x^p^(x^ 

p 2 p"(x^p( 5 ) - 99( 3 ) (X£ M 4 ) (X ? ) ' 

and 

1_ _ ^(x|)^(xg) - y?(xg)y( 3 )(xg) 
p 4 (x^( 3 ) (x ? ) - <p& (xeW(xt) ' 

In particular, p and V are uniquely determined by the point where and its derivative 
vanish. 

Let us now write f( = g + h with g, h positive positive definite (either I? or continuous). 
Then, from Theorem 13.21 g and h have at most 4 complex zeroes and, from the proof of that 
theorem and the remark following it, we know that g = A/e for some f = pe** € C and A > 0. 
But < g = < implies that g must also have a double zero at x^. Hence, by the 
previous argument, we must have £ = £, and g = Xf. 
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4. HERMITE FUNCTIONS 

In this section we restrict our attention to the one-dimensional situation. 

4.1. Preliminaries. 

Hermite functions are functions of the form P(x)e~ Xx with P a polynomial and A > 0. 
They satisfy many enjoyable properties, in particular, they provide the optimum in many 
uncertainty principles, that is, their time-frequency localization is optimal {see e.g. |FSl IHJ| 
and the references therein). 

Let us define the Hermite basis functions by 

nl/4 

h k {x) = ^^ =e ™V e - 2 ™ 2 , k = 0,1,2,... 

It is well-known that (/ifc)fc=o,iv- f° rm an orthonormal basis of L 2 (1R), and hk(x) = CkHj t (x)e~ nx2 
with Hk a real polynomial of degree k with highest order term (2^JlTx) k and a normalization 
constant that is not relevant here. A simple computation shows that 

(8) H (x) = 1, H 2 (x) = 4irx 2 - 1 and H 4 (x) = (4irx 2 ) 2 - 6(4vrx 2 ) + 3. 

Finally, let us recall that the Hermite basis functions are eigenvectors of the Fourier trans- 
form: hk = (—i) k h k . It immediately results that if the Hermite function H(x) = P(x)e~ 7TX 
is positive positive definite, then the degree of P is a multiple of 4. Indeed, H has to be 
even, so that P is real even. Then, if we expand H in the Hermite basis, we get that 

k 

H(x) = a j^2j (x)e~ wx where the ay's are real and ^ 0. The Fourier transform of H 
3=0 

k 

is then given by H{$) = (— 1 ) 3 ' ctj Hyj (£ ) e -7r ^ . But, the highest order terms of the polyno- 

j=0 

mial factor of H and H are then respectively ak(^7TX 2 ) k and (— l) fc afc(4-7r^ 2 ) fc . Checking that 
H and H stay non-negative when x and £ go to infinity suffices to see that k is even. 

4.2. Gaussians are extremals. 

Let us now turn to properties of extremals among Hermite functions. The first result is a 
simple consequence of Hardy's Uncertainty Principle. 

Proposition 4.1. 

Let A > and P be a polynomial of degree N. Assume that f(x) = P(x)e~ XlTX is a positive 
positive definite Hermite function. Then 

1(f) C {Q(x)e~ XlTx2 , Q a polynomial of degree < N}. 

In particular, if f(x) = e~ Xnx2 , then f is an extremal ray generator in 

The second part of this proposition (and its proof) seems well-known, see [Bor] . The proof 
given here is only a slight improvement. 
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Proof. The second part of the proposition immediately follows from the first one. 

Let P be a polynomial of degree N such that f(x) = P(x)e~~ X7TX G Sl + and assume that 
P(x)e~ X7rx2 = g(x) + h(x) with g,h £ 0+. Then, as h(x) > 0, < g(x) < P(x)e" A7rx2 . 
Further there exists a polynomial P of degree N such that P(^)e _7r ^ / A = /(£) = + 
so that, as h > 0, < 5(C) < P^e" 7 ^. 

In particular, there exists a constant C such that < C(l + IxD^e - '* 47 ™ and \g{Cj\ < 

C(l + I^D^e"^ 2 ^. From Hardy's Uncertainty Principle ([Hi] see e.g. [FSl E3] and [BlXTl 
IDell IDe2| for generalizations), there exists a polynomial Q of degree at most ./V such that 
g{x) = Q(x)e~ Xnx2 and therefore h(x) = (P - Q)(x)e~ Xnx2 . □ 

Let us conclude with the following lemma: 
Lemma 4.2. Let A > and P be a polynomial of degree AN > A and assume that f(x) = 
P(x)e~ nXx2 is a positive positive definite Hermite function. Write = i- > (£)e~ 7r £ 2//A . 
- If f is an extremal ray generator, then either P or P has at least 4 real zeroes. 

— If P or P has AN real zeroes, then f is an extremal ray generator. 

Proof. It is enough to prove the lemma with A = 1. Assume that P and P are non-negative 
and have no zeroes, then there exists a constant C > such that, P{x) — C and P{C) — C are 
still non-negative. It follows that p ( x ^~ c e -7 ™ 2 and its Fourier transform, P ^\~ C e - ^ 2 are 
both positive. As f(x) = p ( x ^~ c ' e~ vx2 + p ( x ^ +c e ~ nx2 we conclude that / is not an extremal 
ray generator. It follows that P or P has at least one zero on the real line. As they are 
positive polynomials, any real zero must have even multiplicity. As these polynomials are 
even and non-zero at the origin, we get that P or P indeed has at least 4 real zeroes. 

For the second part, it is enough to consider the case when P has 4A^ real zeroes. If 
f = g + h with g,h € /(/) then, from Proposition 14.11 we get that g(x) = Q(x)e~' KX and 

2 

h(x) = R(x)e~ wx with Q, R polynomials of degree at most AN. As/i>0, < g < f thus 
< Q < P. It follows that a real zero of P is also a real zero of Q. As Q has not higher 
degree than P, Q = cP with < c < 1 and then g = cf, h = (1 — c)f. □ 

4.3. Extremal ray generators among Hermite functions of higher degree. 

2 

From the previous section, we know that if a Hermite function f{x) = P(x)e~ wx is positive 
positive definite then the polynomial has degree a multiple of 4. Moreover, the interval 
generated by / consists of Hermite functions of not higher degree. Thus, if there exists a 
polynomial P of degree Aq such that P{x)e~' nx is positive positive definite, then there exist 
extremals of the same degree. Indeed, we just have to consider the finite dimensional cone 
of positive positive definite Hermite functions of degree Aq, which is then non-empty, and is 
thus the positive span of its extremal rays. 

We will now characterize all positive positive definite Hermite functions of degree 4 and 
the extremals among them. It is enough to consider / of the form f(x) = (Hq(x) + 2aH2(x) + 
bH4(x))e~ 7TX . Using the fact that the Hermite basis consists of eigenvalues of the Fourier 
transform, we get that /(£) = (flo(f) - 2aH 2 (Q + bH 4 (£))e-^ 2 . 

We thus aim at characterizing a, b for which P{x) = Hq(x) + 2aH2(x) + bH^x) > 0. But 
H Q (x) + 2aH 2 (x) + bH 4 (x) = 1 + 2a + 36 - 8vr(a + 3b)x 2 + I6bir 2 x 4 . Setting X = Airx 2 , we 
thus ask whether P(X) := 1 + 2a + 3b - 2(a + 3b)X + bX 2 > for all X > 0. 
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The first condition is that P(0) ^0 that is 1 + 2a + 36 > and that limx^+oo P(X) > 0, 
that is 6 > 0. Next, we want that P has no single root in ]0, +oo). Thus, either a + 36 < 
or (a + 36) 2 — 6(1 + 2o + 36) < which we may rewrite as 

(9) (a + 26) 2 + 2(6-I) 2 <I. 

This is the equation of an ellipse £+, which passes through the point (0, 0) where it is tangent 
to the line of equation 6 = and through the point (—1, 1/3) where it is tangent to the line 
of equation 1 + 2a + 36 = 0. 

In this form, we see that the set T> + of all (a, 6)'s for which P > is the union of the ellipse 
of equation ([9|) and the triangle formed by the lines 6 = 0, a + 36 = and 1 + 2a + 36 = 0, 
that is the triangle with vertices (0,0), (—1,1/3) (—1/2,0). 

If we want / to be positive positive definite, then both (a, 6) and (—a, 6) have to belong to 
T>+. It is easy to see that this reduces to the intersection of the two ellipses. More precisely 
we get : 

Proposition 4.3. The set of all Hermite functions of order 4 that are positive positive definite 
is given by f a ^ := (Hq + 2aH2 + bHi)e~ 1TX where (a, 6) belongs to the region T> parametrized 
by the following : 

- either a > and (a + 2b) 2 + 2(6 - 1/4) 2 < 1/8. 

- or a < and (-a + 26) 2 + 2(6 - 1/4) 2 < 1/8. 

Moreover f a ^ is an extremal ray generator of Q + if and only if (a, 6) € &D, where dT> is 
the boundary ofT>. 



Figure 1. The set of all (a,6)'s for which f a ^ is positive positive definite. 

The last part results directly from the fact that D is convex and that all its boundary 
points are points of curvature (excepted two) and are thus extreme points of T>. The fact 
that fafi is then an extremal of Q + is a direct consequence of the previous discussion. 

It should also be noted that if f a 6 is an extremal ray generator, then 
— either f a ^ has only real zeroes and f a ^ is positive, in which case we will say that f a ^ is 
of the time type 

_ or fafi has only real zeroes and f a ^ is positive, in which case we will say that f a ^ is of 
the frequency type. 

N 

Corollary 4.4. Assume that f ai ,bn ■ ■ ■ , fa N ,b N are a ^ °f ^ e time type, then f = JJ/ ai ,6j is 

i=i 

an extreme ray generator in C. 
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If fa lt b t , ■ ■ ■ ,fa N ,b N are all of the frequency type, then f = * i=1 f ai ,bi is an extreme ray 
generator in C. 

N 

Proof. Let us first assume that all the /a^'s are of time type and consider / = JJ/a 4 ,6i- If 

i=l 

this is extremal, then / = *fLif-ai,bi is also extremal. 

Let us write / = P^ N = g + h with P a polynomial, 7 the standard gaussian, g,h G 
From Proposition 14. 1\ we know that both g and h are Hermite functions, g = G~f N and 
h = Hj N with G, H positive polynomials of degree < AN. 

But, ash>0,0<G<P and every real zero of P is thus a real zero of G. As P is assumed 
to have AN real zeroes, G has 4iV real zeroes. As G is of degree at most AN, G = XP thus 
g = Xf and h = (1 - A)/. 

If all the /a 4 ,f) 4 's are of frequency type, the same argument applies on the Fourier side. □ 



5. Some important classes of positive positive definite functions 

In this section, we will show that several criteria available in the literature are particular 
cases of Choquet representation of the cone Q^. We will appeal several times to lemma [2~T1 

5.1. Non-negative convex functions. 

A well-known class of positive definite functions is that of Poly a type functions (see [Po| ) . 
More precisely, / is said to be of Polya type if / is even, continuous, convex on [0, +00) and 
f{x) — > when x — > 00. It is well-known that / is of Polya type if and only if there exists a 
positive bounded measure v on [0, +00) such that, for all 

(10) f(x)= / (l-\x\/t) + dv(t). 

Jo 

Note that Lebesgue's Theorem ensures that a function / given by (jlOH is continuous as soon 
as the measure v appearing in that formula is bounded. 

Further, one may show that / is a Polya function then / is a Fourier transform of a measure 
of the form p(£) d£ with p continuous on R \ {0}. More precisely, v and p are related by 

/ sin(7rtf ) \ 2 , , . 

Hr 2 ) td " (t) 

where c is some constant (and = 1). 

Note that everything is straightforward when / has sufficient smoothness and decay, using 
integrations by part. For a simple proof in the general case and further references, we refer 
to |Saj . 

Finally, an easy computation shows that (1 — = cx[~t/2,t/2] * X[-t/2,t/2]( x ) which is 

an extremal ray generator according to Theorem 13.21 We may thus rewrite (|10p in the form 

r+00 

(11) f( x )= X[-t/2,t/2]*X[-t/2,t/2](x)du{t). 

Jo 

In this form, we immediately see that we are in the situation of Lemma 12.11 



p(0 
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Proposition 5.1. 

Let f be a Polya type function. Then f is positive positive definite, and it is an extremal ray 
generator if and only if f = X[-r/2,r/2] * X[-r/2,r/2] f or some r > 0. 

Example. 

- ip(t) = ^iqrjjp and ip(t) = n^rms , a > are of Polya type and are thus not extremal ray 
generators. 

- It is easy to show that e~^ P satisfies the hypothesis of the theorem for < p < 1, in 
particular, they are not extremal ray generators. We will show below that this stays true for 
1 < p < 2 as well. 

5.2. Generalizations by Gneiting. 

Polya's criterion has been extended recently by Gneiting [GnTJ IGn2j and by C. Hainzl and 
R. Seiringer |HSj . For instance, [Gnl] characterizes those functions which can be written in 
the form 

r+oo 

(12) VO) = / w(xt)dfi(t) 

Jo 

where w{x) := (1 — x 2 )+ * (1 — x 2 )+ and ^ is a positive bounded measure on [0, +oo). Note 
that uu is the positive positive definite function introduced by Wu in the study of radial basis 
function interpolation and that w(x) = c(l — |x|)^(l + 3/2|x| + x 2 /4) where c is a constant. 
This function also plays an important role in the work of Wendland [We| on positive definite 
functions with compact support of optimal smoothness. 

More precisely, Gneiting showed that a function ifj(t) = (f(\t\) can be written in the form 
(|12p if and only if ip satisfies 

(i) ip is twice continuously differentiable, 

(ii) (/?(0) > and tp(x) — > when x — > +oo, 

(hi) -(Vtip"(Vt)-^'(Vi)) is convex. 

Note that such a function is necessarily non-negative. 

The same proof as in the previous section then shows the following: 

Proposition 5.2. 

A function ip of the form (|12p (or equivalently a function ip given by ip(t) = with (p 

satisfying i) — Hi) above) is positive positive definite, and it is an extremal ray generator 
in the set of positive positive definite functions if and only if there exists c > such that 
tpit) = w(ct). 

Example. 

As already noted by Gneiting, the following functions satisfy (|12p and are therefore not 
extremal ray generators: 

^(*) = -, , urn , < p < 1.877..., ii) ${t) = (1 + 7|t| + t 2 ) exp(-|t|), < 7 < 1/4 
1 + |t| p 

Hi) i/>(t) = (1 - \t\f + (l + 3|t|) iv) i/>(t) = (1 - \t\ x f + 
The first one is known as Linnik's function and is of Polya type for (3 < 1. The third one has 
been introduced by Wendland and has interesting optimal smoothness properties. The last 
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one is called Kuttner's function and seems to be Gneiting's original motivation in the above 
characterization. 

Finally, i/j(t) = exp(— \t\@) satisfies the hypothesis of the theorem for f3 < 1.84170 thus 
improving the domain of non-extremality found in the previous section. 

5.3. Bernstein functions. 

Let us recall that a non-negative function g on [0, +oo) is called completely monotonic if it 
is infinitely differentiable on (0, +oo) and, for all k E N and x £ (0, +oo), 

(-i)V fc) (*)>o. 

Completely monotonic functions have remarkable applications in different branches. For 
instance, they play a role in potential theory [BF], probability theory [Bon} IPe ] IKi] . physics 
|Daj . numerical and asymptotic analysis (KJ IWim] . and combinatorics [BalJ. A detailed 
collection of the most important properties of completely monotonic functions can be found 
in |Widl Chapter IV], and in an abstract setting in |BCR| . 

The celebrated theorem of Bernstein (see |Be|, Chapter III Setion 2], [Fil Chapter 18, 
Section 4] or |WidI page 161]) states that every completely monotonic function which is 
continuous at zero is the Laplace transform of a positive bounded measure on [0, +oo). In 
other words, g is completely monotonic and continuous at if and only if there exists a 
(necessarily unique) positive bounded measure [i on [0, +oo) such that, for every x > 0, 

r+oo 

(13) g(x)= / e- tx dfx(t). 

Jo 

A particular case of Lemma l2.1( of which the first part is well-known, is then the following: 
Proposition 5.3. 

Let g be a non-negative, completely monotonic function, continuous at 0. Let f be defined on 
^ d by f( x ) = 9(\ x \ 2 )- Then f is a positive positive definite function and is an extremal ray 
generator if and only if f is a Gaussian. 

Example. 

— Let A > and < a < 1. It is well-known and easy to show that g defined by g(t) = e~ xt ™ 
is completely monotonic. It follows that e~ x ^ P is positive positive definite for < p < 2 and 
an extremal ray generator only for p = 2. 

- Let a/0 and > 0. It is easy to show that g(t) = (t + a 2 )^ 13 (the so-called inverse 
multiquadrics) is completely monotonic. It follows that (x 2 +a 2 )~ /3 is positive positive definite 
but not an extremal ray generator. 

- Recently, the so-called Dagum family Dp^(x) = 1 — ( \ has been introduced in 

[PMZP, MPN] where it was shown that the Dagum class allows for treating independently 
the fractal dimension and the Hurst effect of the associated weakly stationary Gaussian RF, 
by using the procedure suggested in [GSJ. In [BMP], the authors further showed that for 
certain parameters, this function is completely monotonic. 

Note that completely monotonic functions allow to construct radial positive definite func- 
tions in any dimension. There is a converse to this. More precisely, let /o be a continuous 
function on [0, +oo) and define fd on M d by fd(x) = /o(|ic|). Observe that if fd is positive 
definite then fd' is also positive definite for any d' < d. A famous theorem of Schoenberg [Sc] 
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(see also [SvP] for a more modern proof and further references or |Bax] for another proof) 
states that if fd is positive definite for any d, then fo(x) = g(x 2 ) with g completely monotonic. 

6. Conclusion 

In this paper we have investigated the extremal ray generators of the cone of continuous 
positive positive definite functions for which we have described two important classes of such 
extremals. 

At this stage, we would first like to stress that, up to minor modifications, most of our 
results stay true for positive positive definite L 2 functions or even tempered distributions. 
Let us recall that in this case, one can not define positive definiteness via (I14p . One thus 
replaces this condition by one that is equivalent for continuous functions. To do so, notice 
that if / is continuous, then / is positive definite if and only if, for every smooth function 
$ G Q(R d ) in the Schwartz-class, 

(14) f [ f(x-y)<£>(x)¥&)dxdy>0. 

Note that this makes sense if / is only assumed to be in L°°, but in that case it is well- 
known (e.g. \Be\ \Chl\ IPhj ) that implies that / is almost everywhere equal to a continuous 
function. Further, we may extend (|14j) to T £ 13' and take this to be the definition of a 
positive definite tempered distribution : a tempered distribution T is positive definite if for 
every $ G 8(M d ), (T, $ * <J>*) > where $*(x) = x). Schwartz [Sch| extended Bochner's 
theorem by characterizing positive definite tempered distributions as the Fourier transforms 
of positive tempered measures, that is (T, <&} = J Rd 3>(£) d/i(£) where fi is a positive measure 
such that, for some a > 0, J R d(l + |£ 2 |)~ a d/i(£) < +oo. A proof based on Choquet Theory 
may be found in [To] . As a particular case, one gets that if T £ L 2 , then T is positive definite 
if and only if T > 0. For further extensions of the notion of positive definiteness, we refer to 
[Ste] . 

Unfortunately, several open problems remain, but some features seem to have become 
clear: 

— A full classification of the extremals is probably impossible and no reasonable conjecture 
can be stated at this stage. For instance, the existence of an extremal compactly supported 
function with Fourier transform having complex zeroes leaves little hope for a full description 
of compactly supported extremals. As for the extremals inside the Hermite class, we have 
a full description only up to degree 4. For higher degree, even though we have been able to 
construct extremals, we are far from a full description. The main reason for this, is that we 
are unable to decide the answer of the following question: 

Question 1. 

a) Is the product of two extremals an extremal? 

b) What about the convolution of two extremals (if it makes sense)? 

We have no hint of what the answer might be and will therefore not propose a conjecture. 

— There are many criteria allowing to decide whether a function is positive definite. As we 
have seen, these criteria actually characterize the functions that are mixtures of the scales of 
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a single extremal inside the class of positive positive definite functions. We are thus tempted 
to ask for more criteria, for instance: 

Question 2. 

a) Find a criterion that allows to decide whether a function is the mixture of the scales of 
m a * m a , where a may, or may not be fixed. 

b) Find a criterion that allows to decide whether a function is the mixture of extremals of 
scales of extremal Hermite functions of order 4. 

c) Find a criterion for positive dehniteness that is not given in terms of a mixture of positive 
positive definite functions, thus also characterizing positive definite functions that are not 
necessarily non-negative. 

— Finally, all the extremals we found, except the gaussians, either have at least a zero, or 
their Fourier transform has at least a zero. Borisov |Bor| conjectured that the only extremals 
m L 2 (M)nC(l) that are strictly positive and that have strictly positive Fourier transform 

2 

are the gaussian functions e~ a7TX , a > 0. We are inclined to believe that this is false. This is 
linked to Question 1 above which leads us to the following more precise question: 

Question 3. 

2 

Consider the two functions f = X[-i/2,i/2] * X[~i/2,i/2] an d 7 = e~ nx , and recall that they 
are extremal positive positive definite functions in Qf. Is f * 7 extremal in Qf and, if so, is 
lif * 7) extremal in Q^? 

A positive answer to the second question would of course provide a counterexample to 
Borisov's conjecture. 

Overall, we hope this paper will help reviving P. Levy's request to study the cone of positive 
positive definite functions. 
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